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Abstract 
Enrico Fermi estimated the yield of the Trinity test to be about 10 kilotons by dropping small 
pieces of paper and observing their motion in the blast wave. This is about half the radiochemically 
derived value of approximately 21 kilotons that necessarily includes thermal and nuclear radiation 
that do not contribute to the blast. Although this story is classic, there appears to be no account 
of how he related his observation to the yield. This note attempts to reconstruct how he might 


have done so. 


I. INTRODUCTION 


The official history of Los Alamos! and many popular accounts describe how Fermi, at 
Trinity, the first nuclear test, conducted July 16, 1945, estimated its explosive yield with a 
simple experiment. He dropped some small scraps of paper before the blast wave passed by 
and observed how far they were displaced. From this, he was able to estimate its explosive 
yield as 10 kt. This is about 40% of the modern estimate”, but some portion of that was 
radiated and did not contribute to the blast wave. Fermi did remarkably well with a very 
rudimentary experiment, but neither he nor subsequent writers appear to have explained 
how he inferred the yield. 

The light scraps of paper that Fermi dropped moved with the air because of their low 
ballistic coefficients (the ratio of mass to the product of drag coefficient and area projected 
perpendicular to the velocity through the fluid); they traced the motion of the air in the 
blast wave. The text of Fermi’s memorandum may be found on a number of web sites, and 
it is reproduced in Fig. 1. This note attempts to reconstruct his reasoning. 

Fermi measured the positive phase air displacement D. In the limit of a weak blast wave, 
D is directly related to the positive phase impulse per unit area Z, and measuring D may 


be the best way of measuring Z: 
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where pco is the acoustic impedance of the air. 


What was the understanding of blast wave physics at the time of Trinity? The theory was 
developed in some detail in the 1944 Los Alamos document LA-165", and fitted to empirical 
data on small (mostly 67 Ib) charges in LA-316° (June 25, 1945, three weeks before Trinity). 
Fermi’s responsibilities as a nuclear physicist and experimentalist were elsewhere, but he was 
also a leading theorist and likely familiar with these results. However, the results of? do not 
extend to overpressures below 1 lb/in?, approximately ten times those at Fermi’s observation 
point, and its results for impulse (directly related to the displacement Fermi measured) for 
yields ~ 10,000 tons relevant to Trinity do not extend to ranges > 2 km, in contrast to the 
16 km range to Fermi’s observation, and its results for low impulses are obtained for low 
yields but at the unscaled 100 foot height of burst. Scaling from the May 7, 1945 100 ton 


high explosive calibration and rehearsal shot® might have been possible, but no quantitative 
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pieces of paper before, during and after the passage of the blast 
Sineg at the time, there was no wind I could observe very distinctly 
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On the morning of the 16th of July, I was stationed at the Base 
Camp at Trinity in a position about ten miles from the site of the explosion. 

The explosion took place at about 5:30 A.M. I had my face protected 
by a large board in which a piece of dark welding glass had been inserted, 


"My first impression of the explosion was the very intense flash of light, and 


a sensation of heat on the parts of my body that were exposed. Although I did 
not look directly towards the object, I had the impression that suddenly the 
countryside became brighter than in full daylight. I subsequently looked in 
the direction of the explosion through the dark glass and could see something 
that looked like a conglomeration of flames that promptly started rising. 
After a few seconds the rising flames lost their brightness and appeared as 
a huge pillar of smoke with an expanded head like a gigantic mushroom that 
rose rapidly beyond the clouds probably to a height of the order of 30,000 
feet. After reaching its full height, the smoke stayed stationary for a 
while before the wind started dispersing. it. 

About 40 seconds after the explosion the air blast reached 
I tried to estimate its strength by dropping from about six feet 


ubupey 'W 4 
GlisISSVIONN 
TWNIa 


NOIL 


actually measure the displacement of the pieces of paper that were in 


» process of falling while the blast was passing. The shift was about 23 meters, 


which, at the time, I estimated to correspond to the blast that would be pro- 
duced by ten thousand tons of T.N.T. 
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FIG. 1. Fermi’s original memorandum’. 


data from it are documented inf; its purpose appears to have been to test experimental 
procedure. In addition, the differing specific energies of nuclear and conventional explosives 
would have made scaling uncertain. 

LA-2000", dated August 1947, two years after Trinity, presents a comprehensive review 
edited, and to a substantial extent written, by Hans Bethe, extending the earlier theory and 
presents the results of numerical calculations apparently performed between 1944 and 1947." 
contains not only sophisticated analytic theory, as befits its authors, but also in Chap. 6 
the results of numerical calculations (“IBM solution, results or run”) carried out on the 
primitive computing machines of the time’. Some of those results are used here. 

Most of the energy of an airburst is divided into three parts: radiation, “wasted” thermal 
energy in heated air (partly in the initial fireball and partly in air irreversibly heated by 
the blast wave) and an outgoing blast wave. A small fraction of the energy appears as 
radioactive decay and neutron capture gamma-rays. The displacement of the air (and pieces 
of paper that act as Lagrangian tracers of its motion) far from the explosion includes a 
net outward displacement that accommodates the increased volume of a bubble of hot gas 
(whose subsequent buoyant rise creates the famous “mushroom cloud”), and the oscillatory 


motion of the outgoing blast wave. 


I. THE BUBBLE 


After pressure equilibrium is achieved, the injection of an energy Y, into the air increases 


its volume by 
wy 1) 
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where y is the adiabatic exponent of air and P its initial pressure. At temperatures of 
less than a few thousand K the vibrational modes of nitrogen and oxygen molecules are not 
significantly excited and y = 7/5. At higher temperatures y is less; full classical excitation 
of the vibrational modes reduces it to 4/3, and dissociation and endothermic creation of 
nitric oxide reduce it still further. Most of the volume expansion is attributable to air at 
lower temperatures, so here we adopt y = 7/5. 

For an airburst at the top of a 30 m tower! and an observer at a distance >> 30 m, it 


is possible to approximate the geometry as hemispherical. The effective radius of the hot 
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bubble, defined as the radius of a hemisphere of volume AV, 
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where at the altitude of Trinity of 1.5 km P) ~ 8.5 x 10° dyne/cm?. An observer at distance 
r would see, after all the dynamic motions have settled down but the hot bubble has not 
yet begun to rise, an outward displacement 


= dmr? P) mr? 3r? 
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For Fermi’s observation at r = 10 miles (16 km) 


Yo 
Arx 
r (a) x 10cm, (5) 


where 1 kt (kiloton) = 4.19 x 10"? ergs. For modern estimates” that Trinity had a yield of 


25 + 2 kt this gives 25 cm, which is much smaller than the 2 1/2 meter displacement he 
observed, even were there no radiation losses. The explanation of Fermi’s observation, and 


the reconstruction of his argument, must be sought elsewhere. 


Ill. THE BLAST WAVE 


The sudden creation of a bubble of hot gas produces an outgoing blast wave. This blast 
wave is initially very strong, with overpressure AP = P— P) > Po, where P is the ambient 
pressure and P the pressure of the shocked and subsequently rarefied air. In the strong shock 
and point source limit it is described by an analytic theory developed by von Neumann’, by 
G. I. Taylor! and by Sedov, and described in this issue by Baty and Ramsey'!. This is not 
applicable to the weak shock regime in which Fermi made his measurement. 

LA-20007 also calculates weaker blast waves. However, one crucial element, the fraction 
of the explosive yield that appears in the weak blast wave (|AP| < P) at large distances, 
where Fermi measured it, can only be calculated numerically because it involves the coupling 
of energy through an intermediate strength regime in which no analytic theory is valid. This 
calculation is described in Chapter 6 of LA-2000", but it is difficult to reconcile the results 
there (Fig. 6) with Fermi’s measurement and the known yield? of Trinity. 

Once the blast wave becomes weak, it is, to a good approximation, a sound wave prop- 


agating at the sound speed. The over- and under-pressure AP and the fluid speed decay 
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nearly proportionally to 1/r in spherical or hemispherical geometry. The width of the blast 
wave is nearly constant because the sound speed is nearly its ambient value throughout. 
Idealized displacement and overpressure profiles are shown in Fig. 2. Unlike the bubble dis- 
placement, in the acoustic limit the blast wave displacement returns to zero after the wave 
passes. 

The dependence of shock overpressure on distance from the explosion is shown in Fig. 3 for 
a spherically symmetric explosion into air at ambient conditions (P) = 1.01 x 10° dyne/cm? 
and T = 300K). This Figure combines the numerical (“IBM”) results with the “semi- 
acoustic” and acoustic theory for weak and very weak (acoustic) blast waves. §6.3 of LA- 
2000 states that the calculation was intended to describe a 10 kt explosion but actually 
described a 13 kt (not including thermal or nuclear radiation) explosion. 

For r > Ry the leading edge of the blast wave becomes a weak shock (AP < Po) that 
performs very little irreversible work (x (AP)?). The blast wave propagates at nearly the 
speed of sound, co = 340m/s in air at 15°C, a typical temperature at dawn in July at 
Trinity (LA-2000 takes co = 347 m/s, appropriate to air at 27°C). The passage of a weak 
shock produces a velocity discontinuity 
JEP 


Av ; 
PCo 


(6) 


defining the acoustic impedance pco, where p is the ambient density (very close to the shocked 
density). The pressure profile of the weak blast wave at one radius is shown in Fig. 4, and 
the duration of the positive phase (P > Po) is shown in Fig. 5, both for a yield of 13 kt in 
spherical symmetry. 

In this weak blast wave regime the rate of dissipation of blast energy by conversion to 
thermal “waste energy”, so-called because it does not contribute to the destructive effect 
of the blast, rapidly decreases with increasing r. However, it decreases slowly enough that 
the energy in the blast wave does not approach a constant asymptote, but rather declines 
logarithmically with overpressure and distance (Eq. 429 of LA-165*). 


The energy of a weak hemispherical blast wave whose positive phase has thickness h 


1 Av)? 
Bous = 207?h x 2x 2x 2 x Pl L ; (7) 


One factor of 2 results from the contribution of the compression energy, equal to the kinetic 


energy in an acoustic wave, and another factor of 2 from the contribution of the negative 
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FIG. 2. Idealized blast overpressure and displacement far from the explosion as a function of time 
at one point or of radius at one instant (if plotted as a function of radius the blast is propagating 
to the left and the location of the explosion asymptotically distant to the right). This profile is the 


classic “N-wave” 1? 


, named for its shape; it is produced when an idealized membrane between two 
regions of infinitesimally differing pressures is abruptly removed. This pressure profile is shown 
in Fig. 5.1 of LA-2000. The displacement profile is added here, and both were familiar to Fermi. 
The rarefaction phase of the blast produced by a strong explosion is quite different because there 
are no rarefaction shocks of finite amplitude; Fig. 4 shows the results of an early calculation from 


LA-2000. Far from the explosion the blast propagates as an acoustic wave with nearly constant 


profile, but with amplitude decreasing in proportion to 1/r. 


phase, nominally equal to that of the positive phase (Fig. 2, although Fig. 4 suggests a factor 
of 1.5 may be more accurate), and the factor of 1/3 comes from averaging v? with v varying 
linearly from 0 to Av through the blast wave (Fig. 2). The thickness h is comparable to 


the radius Ry where the shock transitions from strong to weak, scaling œ Q13, but grows 


T 
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OVERPRESSURE IN THE SHOCK IN BARS 


10° 
SHOCK RADIUS IN METERS 


Fig. 6.1 Shock pressure vs distance for nuclear explosion. Time is given in seconds; T = 0 is 
start of IBM run, 0.012 sec after explosion. 


FIG. 3. Calculated shock overpressure vs. radius for a 13 kt explosion in ambient air in spherical 
geometry, corresponding to a 6.5 kt explosion in hemispherical geometry. Note the separate ordi- 
nate labels, arrival times T and Mach numbers M for each curve. At Fermi’s distance of 16 km, 
in hemispherical geometry the overpressure would have been 0.008 bar for a 6.5 kt explosion and 
0.010 bar for 13 kt. AP scales approximately œ Q!/°, where Q is the yield, because the energy 


density is x (AP)? in an acoustic wave and the blast wave thickness scales x Ry x Q!/3. (Fig. 6.1 


of LA-2000) 


logarithmically (LA-165* and Fig. 5) as the weak blast wave propagates. 


Fermi had to estimate h = coAt theoretically because it would not have been easy to 
measure At without photographic or electronic equipment. At AP = 0.008P (16 km 
from a 13 kt explosion in spherical geometry; Fig. 3) the duration of the positive phase 
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Fig. 6.13 Pressure vs distance at a fixed time 


FIG. 4. Calculated blast wave pressure profile far from the explosion for a yield of 13 kt in 
spherical geometry or 6.5 kt in hemispherical geometry. The width and profile P(R) of the blast 
wave vary only very slowly from this radius (the greatest radius for which this is shown in LA-2000) 
outward, while the overpressure AP scales œ 1/r. This profile is not close to the ideal profile of 
Fig. 2, indicating that the model is only approximate. At the shock P/P) abruptly rises from 
unity, but the rise (a vertical line in Fig. 2) is not shown here. According to §6.4 of LA-2000, 
AP = 0.0251P at 6270 m with a positive phase 290 m wide, but the portion of the pressure curve 
with AP > 0.016P) and r > 6200m is not shown. If the shock is actually at 6270 m, as stated, 


the positive phase is 320 m wide. The negative phase is much wider. (Fig. 6.13 of LA-2000) 


At = h/co = 1.14s for the same explosion (Fig. 5). Scaling to 13 kt in hemispherical 
geometry while holding AP constant doubles the effective yield and multiplies Ry and At 
by 21/3. Then At = 1.44s and h = coAt ~ 500m. This is (probably fortuitously) very close 
to the Ry = 490m estimated from Eq. 3. 


At AP = 0.02P), as shown in Fig. 4 at 6 km for 13 kt in spherical geometry, Fig. 5 
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OVERPRESSURE IN SHOCK 


Fig. 6.14 Deviation of the positive phase of a pulse at a fixed distance 


FIG. 5. Duration At of the positive phase as a function of overpressure (original caption has 
“Deviation” for “Duration” ). The blast wave width increases logarithmically even in the acoustic 
regime because the overpressure, although small, is not infinitesimal. The “IBM Run” was the full 
numerical calculation in the moderately strong and strong shock regimes, while the “semi-acoustic” 
theory of LA-2000 calculates the flow to second order in AP/ Pp; a first-order calculation is strictly 
valid only in the acoustic limit. In this limit the positive phase duration At and thickness h are 


related by h = coAt. (Fig. 6.14 of LA-2000) 


indicates At = 1.04s and h = 360m rather than the 320m indicated by Fig. 4 or the 290 
m of §6.4 of LA-2000. This should be taken as a measure of the resolution and precision of 
the calculations. 

From the results reported in LA-2000 (or the more approximate elementary analytic 
result Eq. 3) Fermi had a good idea of the value of h at his observation point for yields in 


the expected range. Fortunately, At and h depend only logarithmically on AP (LA-1654 
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and Fig. 5). In principle, he could have used the computed dependence of h on yield to 
solve interatively for the yield, but this dependence is so weak that such a procedure was 
unnecessary. Measuring a displacement D = 2.5m provided the peak velocity Av for the 
linear pressure and velocity profile of Fig. 2. Then 


= 2Dco 


A 
re a 


= 3.4 m/s. (8) 


This is fairly close to the v = 2.9m/s calculated for AP = 0.010bar (scaling the re- 
sults of Fig. 3 to hemispherical geometry by multiplying AP by 21⁄3) from Eq. 6, taking 
p = 1.0kg/m? at the 1500 m altitude of Trinity and noting that AP scales « Epiast/r?, 
independent of the air density. Substituting h and Av in Eq. 7 yields 


EBlast m 1.5 kt. (9) 


In order to relate his observation to the explosive yield Fermi had to know what fraction 
of the yield appears as blast wave energy at an overpressure of about 0.01 bar. If he had 
relied on the IBM calculation reported in Chap. 6 of LA-2000, shown in Fig. 6, he would 
have been seriously misled. A small extrapolation from these results to AP = 0.01 bar 
would have led to Epiast © 0.016Q, implying an impossible Q ~ 90 kt, far from both Fermi’s 
estimate and the modern radiochemical value of 25 + 2 kt? for the total yield, that includes 
thermal and nuclear radiation. 

Two factors may contribute to this discrepancy. The effects of numerical viscosity were 
likely not appreciated when the IBM calculations were done, and the decline of Egiast/Q with 
increasing range in Fig. 6 may be mostly a numerical artefact. In addition, the “wasted” 
energy contributes to the pressure driving the blast wave, so it is not really wasted. These 
hypotheses could be readily tested with a modern calculation. 

How did Fermi get it right? He could have readily estimated Eglast from Eqs. 7 and 8 or 
equivalent forms. He would have needed a value for h, but that could have been taken from 
the numerical calculations (were they available and known to him in 1945), estimated in 
advance of Trinity as the shock radius Ry (Eq. 3) where AP = O(P), or inferred at Trinity 
from eyeball observation of the duration of the positive phase using his paper Lagrangian 
tracers. He didn’t tell us what he did. 

Fermi would also have needed a value for Epiast/@Q, an essential parameter that the only 


possibly available numerical calculation got wrong. He was likely familiar with the Bethe- 
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Fig. 6.15 The positive impulse and the energy in the blast. IBM run scaled to 40,000 tons in 
free air. (Q is the total energy.) 


FIG. 6. Impulse of the positive phase of the blast and the fraction of the bomb energy Q 
that remains in the blast wave. The impulse is scaled to a 40 kt explosion, but the function 
EBiast /Q(AP/P) relating two dimensionless parameters is approximately independent of the size 
of the explosion (this is only approximate because scaling breaks down for the radiated fraction of 


Q). (Fig. 6.15 of LA-2000) 


Kirkwood-Penney result (Eq. 429 of LA-165*) 


EBlast £5 1 ae 
Q 3,/In P/AP 
which, combined with his empirical value for Egst (Eq. 9) immediately leads to his yield 


estimate Q ~ 10 kt. 


0.15, (10) 


If the IBM calculations had been done by the time of Trinity, either Fermi was unaware 


of them or wisely chose to use the analytic theory rather than those apparently quantitative 


12 


computational results. 


IV. CONCLUSION 


At the time Fermi wrote his famous memorandum, the result (Eq. 10) of LA-1654 was 
likely well-known at Los Alamos, so that no explanation was necessary. This theory is 
no longer common knowledge, even though the original documents are readily available to 
anyone. 

Fermi estimated the yield of Trinity to be 10 kt, about 40% of the modern value, from 
the motion of Lagrangian tracers, his famous scraps of paper. Significant additional energy 
was radiated, but how much depends on the physics and chemistry of hot air (in particular, 
the thermodynamics and opacity of oxides of nitrogen) and requires numerical methods to 
calculate. Fermi’s result was therefore a lower bound, and was sufficient to declare Trinity 


a SUCCESS. 
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